Effect of Ridge- Ridge Interactions in Crumpled Thin Sheets 
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We study whether and how the energy scalings based on the single-ridge approximation are revised 
in an actual crumpled sheet; namely, in the presence of ridge-ridge interactions. Molecular Dynamics 
Simulation is employed for this purpose. In order to improve the data quality, modifications are 
introduced to the common protocol. As crumpling proceeds, we find that the average storing energy 
changes from being proportional to one-third of the ridge length to a linear relation, while the 
ratio of bending and stretching energies decreases from 5 to 2. The discrepancy between previous 
simulations and experiments on the material-dependence for the power- law exponent is resolved. 
We further determine the averaged ridge length to scale linearly with the crumpled ball size R, the 
ridge number as 1/i? 2 , and the average storing energy per unit ridge length as iy/£ 2 - 364 ~ 2 - 487 . These 
results are consistent with the mean- field predictions. Finally, we extend the existent simulations 
to the high-pressure region for completeness, and verify the existence of a new scaling relation that 
is more general than the familiar power law at covering the whole density range. 

PACS numbers: 62.20.F-, 46.32. +x, 89.75.Fb 



Although crumpling is ubiquitous and simple to enact, 
intense researches to understand its complexities only 
began in the last thirty years or so. It has not only 
become relevant to cutting-edge technologies like utiliz- 
ing crumpled graphene sheets [1 to harvest energies by 
converting motion into electricity, but is of interest to 
the general phenomenon of condensation of many out- 
standing problems in physics [2]. However, some prop- 
erties of crumpling remain unresolved. For instance, as 
the developable cones [3 -6 increase in number and form 
the familiar network [7HS] of ridges and vertices, scien- 
tists still do not know how they collaborate to produce 
stunningly simple power laws between the crumpled ball 
size R and crumpling force F pJ3UT9] . R ~ F~ a and for 
the occurring frequency[20-22 of different noise inten- 
sities. One challenge is to understand how the crum- 
pled sheet constructs spontaneously a highly porous and 
yet robust structure [23- 25 . Proper theoretical tools are 
also in demand to tackle the complex many-body in- 
teractions among these deformations, especially the self- 
avoidance that plays an important role [14] at creating the 
glass-like interior after a series of highly non-equilibrium 
processes [26 , similar to the random packing in a golf ball 
basket or salt jar. 

In view of the theoretic inability, Molecular Dynamics 
Simulation becomes a powerful tool to obtain informa- 
tion such as the three-dimensional distribution of ridges 
and facets before unfolding, energy storage in each ridge, 
and how the increase of crumpled ball density affects 
these quantities. We are interested in the effect of ridge- 
ridge interactions on the energy scaling predicted by the 
single-ridge approximation [2j [7J [9]. Attentions are also 
paid to resolve major disagreements between previous 
simulations [14, 18J and experiments p~6| [T9 ] . In the mean 
time, we push the simulation to a more-time-consuming 
region of large densities to verify the existence of a new 



scaling law[T9|. Finally, we study how the ridge number, 
averaged ridge length, and energy cost per unit length 
evolve with the crumpled ball density, and compare these 
relations with the mean-field predictions. 

Our simulations follow a protocol which models a 
thin sheet by a triangular lattice with bond length r = 
1[HJ [H HE]. Bending and stretching moduli[12 , k B 
and ks, are imposed on a two-dimensional circular mem- 
brane with radius Ro=130 lattice sites, which composes 
of iVbead = 62143 beads. Simulation on a larger sys- 
tem of Ro=500 sites, for test, has been checked not to 
alter our conclusions. The WCA[27] Lennard- Jones po- 
tential is used to insure that no bead can penetrate each 
other. Crumpling force is simulated by a collapsing im- 
penetrable spherical wall that clothes the membrane. In 
order to compare with the experiments, a realistic value 
of 1/3 is assigned to the Poisson ratio and the relation 
kB/ks — 3h 2 /32 is used to constrain ks with each choice 
oiks- The bead diameter is held constant, h = 0.9. We 
also include plasticity by halving the magnitude of ks 
beyond a yield angle of 10° and requiring the strain to 
relax with the original k B - 

We first verify that our simulation reproduces the cor- 
rect mechanical response. The data plotted in Fig. [I] 
show how the (dimensionless) density of the crumpling 
sheet L>i = 7V bead [47r(/i/2) 3 /3]/[47r^ all /3] varies with 
the (dimensionless) external pressure P* = (P/Y)(Ro/h) 
where R wa i\ is the radius of the collapsing wall and Y is 
the Young's modulus. The reason why we show D\ vs. 
P* here, instead of R vs. P (or F), is that the presenta- 
tion collapses the simulation data obtained from different 
ks onto a master curve, irrespective of whether plasticity 
is included. The results reveal a general scaling relation 
across the whole range of P, which can be reduced to a 
power law Di ~ p*3a/(i+2a) at sma n density. We find 
that the exponent a is 0.25, in good agreement with sim- 
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FIG. 1: Dimensionless density D\ vs. dimensionless pressure 
P* for circular elastic and plastoelastic sheets of different ks- 
Experimental data 19 for Al and HDPE have been superim- 
posed on the plot for comparison. The power-law region is 
marked on the left-hand side of the vertical dashed line. 

illations obtained by other groups [21 [18]. However, dis- 
crepancy is observed when we superimpose experimental 
data [19] on the same figure for comparison. The den- 
sities for Al and HDPE are significantly larger than the 
simulation ones. 

The reason leading to the discrepancy is mainly due 
to the coarse-grained triangular lattice model used here. 
The actual volume of sheet is underestimated in the nu- 
merator of Di, which should be ttRqJi where h is the 
thickness of sheet. Also, the volume of the crumpled 
sheet in the denominator of Di is overestimated. A 
correct value can be calculated from simulations which 
is smaller than the volume of spherical boundary V = 
47ri?^ all /3 by the amount of void AKoid (the black vol- 
ume in Fig. [2^b)) . However, since a real paper can not 
cut into itself, the beads that are increasingly wedged in 
the interstitial holes of the triangular lattices as crum- 
pling progresses, as sketched in Fig. |2|a), should not be 
allowed. Had this constraint been imposed, the crumpled 
ball size would swell by approximately the total amount 
of the wedged volume AV we d ge d' 

Therefore, an appropriate dimensionless density should 
be denned as D 2 = (7r^/i)/(V-AKoid+AV W ed g ed). Fig- 
ure [2|c) presents the plot D 2 vs. P* for the different sets 
of simulation data together with the experimental ones, 
where the ratios for the volume differences, AV VO id/V 
and AV we d g ed/V are also shown in the inset for illustra- 
tion. We observe that the sets of simulation data are 
now collapsed with the experimental ones, while retain- 
ing the scaling relation. In the mean time, the power-law 
region shrinks and the exponent reduces to around 0.22. 
The latter is not hard to understand because AV vo [a/V 
in the inset of Figj2|c) dominates at small density and 
raises D\, while AVvj e & ge &/V rules at large density to re- 
duce D\. Both effects combine to level off the line that 
represents the power law in Figjl] 

Please notice that experimentalists usually focus on 
the response in a relatively large density region where the 
power-law relation is not so well-established, as shown in 
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FIG. 2: (a) Schematic plot for a (red) bead partially wedged 
in the interior (the grey prism) of a triangular lattice of 
(blue) beads, (b) Void (in black color) of a crumpled ball 
(c) Dimensionless density D2 vs. dimensionless pressure P* . 
Data are replotted from Fig^ The power-law region shrinks 
from P* < 1.8 to 0.28, and a decreases from 0.25 to about 
0.22. Inset shows that AV vo id/V and AVbe&d/V dominate the 
small and large pressure region, respectively, for the case of a 
ks = 1000 elastic sheet. 

Fig.[2]for Al and HDPE. It explains why a obtained from 
experiments [16\ H9] is usually smaller and non-universal. 
This is because there exists a lower limit to apply exter- 
nal pressure P in crumpling experiments [T6] [T9] where 
high-pressure nitrogen gas was used to provide the am- 
bient pressure. Since the pressure is balanced, the sheet 
needs to be pre-crumpled and packed by a PVC wrap 
with the wrap interior maintained at atmospheric pres- 
sure via a PE tube connecting to the outside of the pres- 
sure chamber. This procedure sets a lower value to P 
and, therefore, it is not easy for experimentalists to in- 
vestigate crumpling in the low-pressure region. 

It is known that the single-ridge approximation [2j [71 
[9] predicts a scaling law for the energy E to create a 
folding ridge on a sheet as a one-third power of the ridge 
length t. If E is partitioned into the bending energy 
Eb and the stretching energy Es, the ratio of the two 
energies is shown to be Eb/ Es = 5. Although the former 
prediction has been checked by simulations for a solitary 
ridge[18 , the two predictions have not yet been verified 
in a realistic case of crumpling where multiple ridges are 
formed and interact between each other. Since the ridges 
are not formed in an equilibrium state, a full relaxation 
of stress is not expected, particularly when the sheet is 
crumpled into a very tight space. 

We apply watershed algorithm to determine the fold- 
ing ridges of sheet. The data are improved manually 
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FIG. 3: (a) Network of ridges on a crumpled elastic sheet 
of k B = 1000 at D 2 = 0.339. (b) Same as (a) but on a 
plastic sheet, (c) N(£) as a function of £. The distribu- 
tion is fitted by N(£) = A/(£y/B) exp{-[ln(£/£ )/ B} 2 } with 
(A,B,£ ) = (429.6,0.66,20.5) for elastic sheets (circles) and 
(449.9,0.48, 19.0) for plastic sheets (squares). 
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FIG. 4: (a) Averaged storing energy E(£) in a ridge as a 
function of £ at D 2 = 0.339 and ks — 1000. Least square fit 
yields 7 = 0.802 and 1.029 for the elastic and plastic sheets, 
(b) 7 vs. 1/D 2 (c) Eb/Es vs. l/£>2 



by connecting broken segments of the ridges obtained 
(see Figj3ja),(b) for example). We observe that the 
distribution of ridge length follows a log-normal dis- 
tribution as shown in Figj3j which agrees with theo- 
retical predictions [2j [7], experiments [28], [29] and other 
simulations [14]. Moreover, we calculate the energy asso- 
ciated to each ridged generated. The mean energy E(i) 
exhibits a power-law dependence P on the ridge length 
£, as shown in FigJSJa). The exponent 7 is then plotted 
in Figj4jb) as a function of density D2. We find that 
the predicted value 1/3 of Witten's single-ridge approx- 
imation is recovered in our simulations at low densities. 
As the density increases, 7 rises and saturates asymptot- 
ically to 1. In the mean while, the ratio Eb/Es drops 
from the single-ridge value 5 to about 2 (see FigQc)). 



Therefore, the weight of energy gradually shifts from Eb 
to Es as D2 increases, whereas both energies are pro- 
moted by the ridge-ridge interactions. Since the depen- 
dence of Eb/Es on D2 cannot be predicted by the single- 
ridge approximation, our results show that the ridge- 
ridge interactions play an important role on the energy 
storage of sheet. 




mean- field picture. Figure [5 
or equivalently £q ~ R. ~" 



FIG. 5: (a) N£ 2 vs. D 2 for an elastic sheet with k B = 1000. 

(b) £0 vs. D2. The power-law exponents are -0.292 and -0.307, 
respectively, for elastic (circles) and plastic (squares) sheets. 

(c) Power law dependence of e on D 2 . The exponents are 
0.788 and 0.829, respectively. 

Having seen how the ridge-ridge interactions revise the 
energy scalings, we now study some properties that turn 
out to obey the prediction of a simple mean-field theory. 
As shown in Fig.[5|a), the ridge number N multiplies the 
square of the mean ridge length £q is essentially constant 
over the whole range of density D 2 . Suppose that the 
ridges are arranged in a square lattice, the area of this 
checkboard is approximately N£q/2, which is consistent 
with the sheet area ttRq = 5.3 x 10 4 . Thus, the fact that 
N£q does not change with D2 can be explained within a 

1 /3 

b) shows that £0 ~ 
lis result can be explained 
simply through a dimensional analysis: £0 scales merely 
with the length scale R of system and not with ks and 
ks, because ridges appear also on a crumpled sheet even 
in the condition ks = ks = 0. 

Furthermore, the average energy per unit ridge length 
follows a power law e ~ £)0.788~o.829 - n ^ low-density 
region D2 < 0.28, as shown in Fig. [5|c), which is 
also the region where R ~ F~ a . It implies that 
the total energy stored in the ridges, N£qe, scales as 

^-1^-3x0.788—3x0.82^ Qr jj-3.364~-3.487. Because of 

energy conservation, the total work W done on the sheet 
during crumpling is equal to the increase of internal en- 
ergy. W can be calculated by F • dR\ and yields a 
scaling #(-V«)+i as R < j r Qj since F ~ R~ l / a . For an 
exponent a ~ 0.22 as obtained in Fig 2jc), the scaling 
reads as W ~ i? -3 54 , which is in good agreement with 
the scaling relation for N£$s here. 

In conclusion, Molecular Dynamics Simulation has 
been employed to study the effect of ridge-ridge inter- 
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actions on a crumpling of self- avoiding sheet. Energy 
scaling based on the single-ridge approximation is valid 
only at low densities. Revision is needed as the ridge- 
ridge interaction intensifies in the high-density region. 
We manage to conciliate the different exponents a be- 
tween simulations and experiments. We investigate how 
the average ridge number, the average ridge length, and 
the average storing energy per unit length vary with the 
density. The relations are found to be consistent with the 



mean-field predictions, in spite of the strong ridge-ridge 
interactions. Particularly, we extend our simulations to 
the high-pressure region and confirm the existence of a 
general scaling relation which is valid for the whole range 
of densities. 
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